Self-Diffusion of Iron in Ll FePd film - as revealed by 
reflectometric methods 



D. G. Merkel, Sz. Sajti, F. Tancziko, M. Major, Cs. Fetzer and L. Bottyan 

KFKI Research Institute for Particle and Nuclear Physics, 
P.O.B 49, H 1525, Budapest, Hungary 



A. Kovacs 

MTA Research Institute for Technical Physics and Materials Science (MFA), 
H 1525, P.O.B. 49 Budapest, Hungary 



A. Riihm, J. Major 

Max Planck Institute for Metals Research 
Heisenbergstr. 3, 70569 Stuttgart, Germany 



R. Riiffer 
European Synchrotron Radiation Facility, BP 220, 38043 Grenoble CEDEX 9, France 



Abstract 

Llo (CuAu(I)-type) ordered FePd as well as FePt and CoPt, due to their high magnetic 
anisotropy, are candidate materials for future ultra-high density magnetic recording. Atomic 
diffusion governs the structural relaxation and associated changes in the physical and magnetic 
properties in these alloys. Such processes involve diffusion paths of a few angstroms. In 
developing modern storage devices it is indispensable to understand these processes. To 
investigate self diffusion, an isotope sensitive technique is needed. Neutron reflectometry (NR) is 
a suitable non destructive method to study self -diffusion in a chemically homogenous isotope 
multilayer with diffusion lengths of the order of a few angstroms. Another recently emerged 
isotope sensitive method for studying self diffusion is Synchrotron Mossbauer Reflectometry 
(SMR) 

Partially ordered and disordered isotope-periodic multilayer (ML) films of composition 
[ nat Fe47Pd53(3 nm)/ 57 Fe47Pds3(2 nm)]io were epitaxially grown on MgO(001)by molecular beam 
epitaxy (MBE). Post grown annealing was performed in UHV chamber at varying temperature 
between 370 °C and 570 °C . The heat treatment time changed from 90 min to 1800 min. 

According to our model, three different micro structural species were identified in the FePd layer 
and were related to the different hyperfine (hf) field distributions in the CEMS spectrum taken 
on the samples (top right figure). The ordered Llo structure attributed to the low hf fild with 
magnetisation pointing along the crystallographic c direction. The intermediate hf distribution 
with rando magnetic orientation represents the disordered fee FePd component and the high hf 
field corresponds to an iron rich environment which is magnetically coupled to the matrix in 
what it is embedded. 

By fitting the total diffusion parameters for all FePd samples obtained from NR and SMR 
evaluation we determined activation energies (Q) and preexponent factors (D°) for the distinct 
environment in the FePd samples. 

We obtained the highest activation energy (1.82+0.38 eV)and lowest preexponential 
factor((5.76±0.35)*10" 14 m 2 /s) for the ordered Ll as oppose to the iron rich regions were the 



activation energy (1.48+0.26 eV) found to be lowest and the preexponential factor the highest 
((1.01±0.61)*10~ 13 m 2 /s). 



Introduction 

Due to the high perpendicular magnetic anisotropy, Llo FePd(Pt), CoPd(Pt) are in great interests 
for future high density magnetic recording media [1, 2, 3], The structure of Llo corresponds to 
CuAu(I)-type structure which includes the alternation of Fe and Pd atomic planes. The formation 
and transformation of Ll phase in such films is controlled by short range (atomic scale) 
diffusion, therefore it is widely studied. In the case of chemical driving force the interlayer 
diffusion can be followed by X-Ray techniques but due to the concentration dependence of the 
diffusion coefficient, it cannot be described by only one Arrhenius equation, even if the 
annealing was performed on increasing temperatures. It is more advantageous if the intermixing 
is followed in only one crystallographic direction in an epitaxial, chemically homogenous alloy, 
where the diffusion doesn't have chemical driving force (besides the isotopic effect) and the 
isotopes of the same elements move according to the same diffusion coefficient. This method 
allows us to measure diffusion on a very short scale. It is very important in the study of 
metastable phases because through the extended diffusion, phase transition occurs hence the 
diffusion parameter changes irreversible. In previous measurements, however the authors don't 
mention, this effect can be observed. In the case of FeZrN the change of conversion electron 
Mossbauer (CEMS) spectrum after heat treatement indicates a phase transformation in the 
sample [4,]. In the case of FePt, according to the paper published by the same group, CEMS 
results show the rotation of the initially out of plane c-axis of FePt into the plane of the sample 
after increasing temperature annealing. According to Rennhofer et al. [5] the diffusion length 
differs in the crystallographic c- and a-b direction the diffusion constant changes not only due to 
the temperature. These systematic errors can only be eliminated (or take into account) if the 
phases or the orientation of the phases are followed by an independent measurement. CEMS is 
an ideal method where Fe containing samples are studied. For this reason in this study we have 
measured CEMS on all samples before and after annealing and the change of the structure was 
taken into account. 



In former studies [4,5] the evaluation of synchrotron Mossbauer spectra (SMR) was performed 
by the normalization of the Bragg peak to the total reflexion peak but Andreeva showed [6] that 
ratio of the intensity of the total reflexion and Bragg peak is very sensitive to the hyperfine 
properties of the sample and also to the conditions of the sample during measurement (0 and tilt 
angle setup respect to the incident beam). For this reason even if the sample is a periodic 
multilayer, it is necessary to know the properties of the sample (namely CEMS spectrum) for the 
determination of the diffusion profile from the SMR spectrum. The error of the peak ratio can be 
extremely large if a capping layer is used or Kiessig oscillation is present due to the buffer layer. 
For the proper evaluation the fitting of the whole curve is required and even better if the promt 
reflectivity spectrum and the delayed spectrum is fitted simultaneously. 

Experimental 

nat Fe47Pd53 7 Fe47Pd53 isotope periodic multilayers with different initial structure were grown on 
x20x2 mm MgO(OOl) substrates by molecular beam epitaxy (MBE) technique. Depending on 
the desired structure the substrate temperature was held at RT or 350°C during growth process. 
For epitaxial deposition 3 nm Cr seed layer was applied on the substrate followed by 15 nm Pd 
buffer layer. Then the bi-layer sequence of 3 nm nat Fe 4 vPd 5 3 and 2 nm 57 Fe47Pd 53 was repeated ten 
times. All nat Fe47Pd 53 and 57 Fe 47 Pd53 layers were prepared by co-evaporation of Fe and Pd at a 
rate of 0.0485 A/s for both 57 Fe and nat Fe and 0.0682 A/s for Pd. To avoid oxidation of Fe a 1 nm 
Pd layer was grown on top of the sample (capping layer). The base pressure of the evaporation 
chamber was 2xl0" 10 mbar which raised 2.8xl0" 9 mbar during the growth process. To achieve 
identical layer structure for the different samples the tooling factors were carefully calibrated by 
test films and the layer thicknesses and composition were determined by x-ray reflectometry and 
Rutherford Backscattering. During evaporation the growing rate was controlled by two 
independent quartz layer thickness monitors. 

-3 

Depending on the experiment the samples were cut into eight (10x5x2 mm for SMR) or four 
(10x10x2 mm for NR) equal pieces. From each set one sample was left in the as deposited state 
and the rest were annealed between 297°C and 530°C for various times between 90 min and 
1800 min. The heat treatment was carried out in UHV conditions. To get the real sample 
temperature, the sample surface was monitored with an infrared heat camera. 



CEMS spectrum was recorded on all samples before and after heat treatment. The experiments 
were performed by using a Co(Rh) single line Mossbauer source with a homemade gas flow 
single-wire proportional counter operating with He gas with 4.7% CH4 extinction gas and at a 
bias voltage of 830+10 V. 

SMR measurements were performed at the ID 18 and ID22 beam line of the European 
Synchrotron Radiation Facility (ESRF, Grenoble, France) in 16 bunch mode. For sequential 
monochromatization of the beam to the 14.4 keV Mossbauer transition ((A,=0. 86025 A) of 57 Fe, a 
Si(lll), then a Si(4. 2. 2)/Si( 12.2.2) double-channel cut monochromator were used. Prompt (non- 
resonant) and Fe delayed (time integrated) nuclear resonant reflectograms were recorded on 
each sample. 

NR curves were recorded at NREX + beamline in (FRM II), Garching. The instrument was used 

o 

in simple reflectivity setup at a wavelength of 4.2A using a supermirror filter for higher 
harmonics suppression. The reflected intensity was detected with a 200x200 mm two- 
dimensional He wire chamber detector. At each reflection angle the counts were summed up in 
a carefully selected region (ROI) to determine the reflected intensity. 

High angle X-ray diffraction was carried out using a D8 Discover type diffractometer (equipped 
with a Gobel-mirror on the primary side) in Bragg-Brentano geometry using Cu Ka radiation 

o 

(A,=1.5415 A). In order to decrease beam divergence, 0.6 mm slits were used after the source and 
in front of the detector. 

Results and discussion 

Fifty samples of layer composition isotope-periodic 57 FePd/ nat FePd multilayer samples have been 
investigated after various retention time and annealing temperature. The interdiffusion of the 
adjacent isotope-periodic layers was determined by fitting the SMR and NR curves (examples of 
which are displayed as inset A, C, and B, D in Fig 3, repectively). In Figure 1 the fitted values of 
the diffusion coefficient are plotted as a function of the reciprocal temperature of annealing 
(Arrhenius plot). 
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Figure 1 Diffuzion coefficients for various annealing temperatures of FePd obtained from NR and SMR 
experiments 

Obtaining activation energy and preexponal factor from fitting the points we get 1.51+0.31 eV 

1 Q 9 1 

for activation energy and and (1.38+7.58)xl0" m s" ) for pre-exponential factor. The points in 
Figure 1 show very diffuse feature. There are orders of magnitude (high above measurement 
error) scattering of the values at a given temperature. This behavior presumes that however the 
samples have the same nominal composition but there can be difference in the inner structure. To 
confirm this, CEMS measurements were taken on all samples. Three distinct microstructural 
species, the low hyperfine (hf) field, a large hf field, and an intermediate hf field species could be 
identified, which were attributed to the ordered Ll , one iron rich nanoncluster phase and the 
disordered fee structural units, respectively. On Figure 2 three selected Mossbauer spectra with 
their hyperfine field (hf) distributions are shown. On the top, the main component is Llo FePd 
(78%) with a 16% disordered fee structure. The presence of iron rich clusters are almost 
neglectable (6%). In the middle the Ll and disordered environments gives almost identical 
contribution (-40%) to the spectrum and the ratio of high hf component is also significant 
(-20%). In the sample which is shown in bottom the high hf field enviroment is completely 



missing, only a little fraction of Ll (13%) FePd and mostly disordered fee components are 
present. 
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Figure 2 Three selected Mossbauer spectra representing higly ordered, partially ordered and disordered 
FePd structures 

These results indicate that the samples should be considered as inhomogeneous. 

If we consider a ID model where the diffusion coefficient (D) is piecewise constant, thus D takes 
an arbitrary positive value D t in volumes of arbitrary size, the diffusion of Fe in the multilayer, is, 
like in a random alloy, a weighed sum of the individual diffusion coefficients of iron in the 
different Fe-environments, , i.e. the total diffusion can be written as follows. 

D = X Ll D LX + Xfcc D fcc + X Fe D Fe 

Where jc L io, #fco*Fe and D LW , D icc , D Fe are the ratios and the effective diffusion parameters of 
the ordered Llo, disordered fee and iron-rich phases, respectively. Expressing the above equation 
we get the following form: 



D(T) = x Ll D° Lln e " +x fcc D° fcc i^ +x Pe D° Pe e'^ 



,l Q L,iQ 



By fitting the diffusion values for all 49 samples obtained from NR and SMR evaluation we 
determined activation energies (Q) and preexponent factors (D°) for the distinct environments in 
the FePd samples. On Figure 3 the measured diffusion and the result of fitting is shown. For 
selected samples the corresponding NR or SMR spectrum is plotted. On the top of the figure the 
composition of the samples are presented. 
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Figure 3 The measured and fitted diffusion coefficients for all FePd samples (bottom). On the top of the 
figure, the relative fraction of the different environments in a given sample is shown. Selected NR and 
SMR spectra after different annealing temperature are shown in the insets 



From fitting the pre-exponential factors and the activation energies founded to be 
D° L1 o=(5.76+0.35)*10" 14 mV 1 , Qlio=1. 82+0.38 eV, D° fcc =(l .32+0.45)* 10" 13 mV 1 

Qf CC =l. 48+0.26 eV, D° Fe =(1.01+0.61)*10" 13 mV 1 , Q Fe =l. 39+0.23 eV for Ll ordered, fee 
disordered FePd and for iron rich regions, respectively. Note, that these values are in the 
crystallographic c direction which is perpendicular to the sample surface. According to these 
results the diffusion in the 57 FePd/ nat FePd is controlled by the disordered fee and the iron rich 
regions. In the ordered Llo regime the atomic movement is almost blocked compared to the other 
structures. 

Summary 

Partially ordered and disordered isotope-periodic multilayer (ML) films of composition 
[ nat Fe 47 Pd 5 3(3 nm)/ 57 Fe 47 Pd53(2 nm)] i were epitaxially grown by molecular beam epitaxy 
(MBE) on Pd(001)(15 nm) buffer and Cr(3nm) seed layer on MgO(OOl). A Pd(lnm) capping 
layer was used to prevent oxidation. Post grown annealing was performed in UHV chamber ("10" 
10 mbar) at varying temperature between 370 °C and 570 °C . The annealing time varied from 90 
min to 1800 min. During annealing the surface temperature of the sample was monitored by an 
infrated thermometer. 

According to our model, three different microstructural species were identified in the FePd layer 
and were related to the different hyperfine (hf) field distributions in the CEMS spectrum taken 
on the samples. The ordered Llo structure attributed to the low hf fild with magnetisation 
pointing along the crystallographic c direction. The intermediate hf distribution with rando 
magnetic orientation represents the disordered fee FePd component and the high hf field 
corresponds to an iron rich environment which is magnetically coupled to the matrix in what it 
is embedded. 

In our model the diffusion (D) piecewise constant, thus D takes an arbitrary positive value Z), in 
volumes of arbitrary size, so the diffusion of Fe in the multilayer, is, like in a random alloy, a 
weighed sums of the individual diffusion of iron of the different Fe-environments. 



By fitting the total diffusion parameters for each FePd samples obtained from NR and SMR 
evaluation we determined activation energies (Q) and preexponent factors (D ) for the distinct 
environment in the FePd samples. 

We obtained the highest activation energy (1.82+0.38 eV)and lowest preexponential 
factor((5.76+0.35)*10" m /s) for the ordered Llo as oppose to the iron rich regions were the 
activation energy (1.48+0.26 eV) found to be lowest and the preexponential factor the highest 
((1.01+0.61)*10" 13 m 2 /s). 
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